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Dynamics of a state of interest coupled to a non-Markovian environment is studied for the first 
time by concatenating the non-Markovian quantum state diffusion (QSD) equation and the Feshbach 
projection operator partitioning technique. An exact one-dimensional stochastic master equation 
is derived as a general tool for controlling an arbitrary component of the system. We show that 
the exact one-dimensional stochastic master equation can be efficiently solved beyond the widely 
adapted second-order master equations. The generality and applicability of this hybrid approach is 
justified and exemplified by several coherence control problems concerning quantum state protection 
against leakage and decoherence. 
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I. INTRODUCTION 

Feshbach projection operator partitioning technique 
(termed PQ partitioning for brevity) allows one to fo- 
cus exclusively on the dynamics of a small subspace in 
the whole Hilbert space of an open or closed quantum 
system [l|, 0. This has been proven to be exception- 
ally useful in dealing with the decoherence suppression 
for quantum states within many interesting contexts in 
physics such as quantum state storage where PQ parti- 
tioning can significantly reduce the resource by focusing 
on a target state rather than on the entire multi-level 
(even infinite level) Hilbert space. The targeted compo- 
nent in the Hilbert space, called P subspace, could be a 
superposition of some or all of the energy levels, while the 
rest of the state space is denoted as Q subspace. Gener- 
ally there is a bidirectional wave-function fiow between P 
and Q subspaces, unless some intervention mechanism is 
implemented. A typical example of this case is the Quan- 
tum Zcno Effect, that may decompose the whole space 
of the system into isolated Zeno subspaces [1, |3| and the 
remaining parts of the whole space. When the system is 
embedded in a dissipative environment, the dissipation of 
the whole system and the leakage of the subspace will mix 
with each other to leave a compelling, yet subtle question 
of the dynamical control of one subspace in open quan- 
tum systems, especially in those with a non-Markovian 
environment. 

A majority of prior research efforts on a non-Markovian 
environment are based on a second-order master equa- 
tion in terms of coupling constant obtained by using 
projection-operator technique Q or the Liouvillean ap- 
proach Q, where the weak coupling regime is typically 
assumed. The work presented in this paper will investi- 
gate the dynamics and control of an arbitrary subspace of 
an open system beyond the weak non-Markovian regime 



by employing the non-Markovian quantum state diffu- 
sion equation [7l-[ll|. which is capable of dealing with 
the strong coupling strength and the arbitrary correla- 
tion function of the environment. The non-Markovian 
QSD equation may be cast into a convolutionless form, 
hence it also serves as a useful tool in deriving the corre- 
sponding exact master equation [l^, HI, [HI . In addition, 
the approximate QSD equation obtained by perturbation 
may also include the contributions from the high-order 
non-Markovian master equation [l3 - [l6j . 



Based on a microscopic quantum dissipative model, 
we will first derive a general one-dimensional dynamical 
equation for a subspace by combining the QSD equation 
and the PQ partitioning technique. The system survival 
probability in the subspace is obtained by the ensemble 
average over the trajectories for the survival amplitude. 
For simplicity but without loss of generality, we consider 
an initial pure state as our subspace of interest (with 
dimension one). 



We show in this paper via the stochastic approach that 
the non-Markovian dynamics of the one dimensional sub- 
space can be effectively controlled by a sequence of pulses 
applied to the system. In particularly, we show that the 
control function c{t) can be chosen as a simple periodic 
rectangular interaction c{t) — ^, where ^I* is the 

interaction intensity, for regions nr— A<t<nT,n^l 
integral; otherwise c{t) = 0, where 1/r is the frequency 



of the pulse and A is the duration time of the pulse in one 
period. The pulses applied here are neither an ideal zero- 
width pulse nor an optimized Bang-Bang control [Hj. 
Yet it is sufficient for our control scheme if the interval r 
is short enough. We find that increasing ^! is also bene- 
ficial to reduce the leakage rate of the subspace. 
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II. FESHBACH PQ PARTITIONING AND 
ONE-DIMENSIONAL STOCHASTIC MASTER 
EQUATION 

We consider an open quantum system with (N + 1) 
normalized base vectors {N is arbitrary) coupled to a 
bath of harmonic oscillators described by the following 
total Hamiltonian (setting h = 1): 



sys 



^{9*xLai+ gxL^'ax) + ^ujxa[ax, (1) 



where H^ys and L are the system Hamiltonian and the 
Lindblad operator, respectively. The system state de- 
scribed by a stochastic wave-function is governed by 
the non-Markovian QSD equation 0, 

idt^t = [Hsys + iLz; - iL^Oit, z*)] = i^cff^t. (2) 

Here is the colored noise arising from coupling to 
the environment such that its statistical mean recovers 
the environment correlation function a{t,s): M[ztz*] = 
a{t,s). Note that 0(t, z*), with its explicit expression 
given below, is the system operator representing the ef- 
fect of the environment. Thus the effective Hamiltonian 
Hcs contains all the information about the open system 
and its interaction with the environment. Whenever the 
operator O is exactly constructed, then effective Hamil- 
tonian is exact. That is, it is directly derived from the 
total Hamiltonian without using any approximations, in 
particular, without Born-Markov approximation. 

The PQ partitioning technique can divide the system 
wave- function (with (7V-|-l)-dimension) into two parts: 
a scalar function Pit) associated with a chosen vector 
denoted by |0) and an A^-dimensional vector Q{t). Note 
that |0) can be an arbitrary component of the system and 
doest not necessarily denote the ground state. With this 
partition, the state and the effective Hamiltonian may be 
written as 



h 


R 


w 


D 



(3) 



Here the 1 x 1-matrix h and the N x iV-matrix D corre- 
spond to the self-Hamiltonians living in the P subspace 
and the Q subspace, respectively. For the effective Hamil- 
tonian, the resulting W and R are not mutually conju- 
gate to each other in general. Consequently, the QSD 
equation ([2]) may be decomposed into two parts: 

idtP = hP + RQ, (4) 

idtQ = WP + DQ. (5) 

The solution for Eq. ([5]) could be formally expressed by 

Q{t) = -i j dsG{t, s)W{s)P{s) + G{t, O)Q(O), (6) 
Jo 

where dtG{t,s) = -iD{t)G{t,s) and G{t,t) = 1. Note 
that the propagator may be written as 



G{t, s) =T^ < exp 



D{s')ds' 



(7) 



where 71- is the time-ordering operator. Substituting 
Eq. ([6]) into Eq. (jlj, we obtain a closed one-dimensional 
master equation for P{t): 



idtP{t) = h{t)P{t) - i / dsG{t, s)P{s) 
Jo 

+ R{t)G{t,0)Q{0), 
G{t,s) = R{t)G{t,s)W{s). 

In case that matrix D is diagonalizable or G{t, s) 



(8) 



E 



N 



n){n\ is a good approximation. 



Eq. can be readily calculated. Especially when 
G{t, s) = 0, a formal solution can be given. 



P{t) 



P{0)~i / ds'R{s')Q{s')e'^o ^'''''■'^ 



'ifo d,sh{s) ^ 

(9) 

Notably, the off-diagonal term R{t) between P and Q 
plays a vital role in the dynamics of P{t) while the local 
control term h{t) only provides a phase factor. 

Equation ([5]) will be the main result in the concate- 
nation of the PQ partitioning and the non-Markovian 
QSD equation described by iJoff- As shown below, 
the one-dimensional master equation can yield some 
very interesting analytical results that are hard to ob- 
tain without invoking certain approximations such as 
weak-coupling approximation. We will now discuss three 
examples, where 0{t,z*) is exact, hence H^e is exact 
and time-local, to illustrate both the free and controlled 
dynamics of states of interest coupled to a zero tem- 
perature heat bath, in terms of their fidelity with re- 
spect to the initial state. For an arbitrary pure initial 
state ipo, it is easily to show that the fidelity J^{t) — 
M[\P{0)*P{t) + gt(0)g(i)|2]. A simpler expression for 
the fidelity is obtained J"(i) = M[|P(t)|2] whenQ(O) = 0, 
which is equivalent to the population or survival proba- 
bility over ipo- Remarkably, in many physically interest- 
ing examples as shown below, the PQ partitioning com- 
bined with the QSD equation can yield simple analytical 
solutions of the fidelity J^{t) . 



A. Dissipative two-level atom model 

The model is represented by Hgys — Eo{t)\0){0\ + 
Ei{t)\l){l\ and L = |0)(1|. The effective Hamiltonian 
in the basis |0), |1) becomes 



HcS — 



Eo{t) 







Ei{t)-iF{t) 



(10) 



where F{t) is the coefficient function in the operator 
0{t,z*) = F{t)L with F{t) = dsa{t,s)f{t,s) and the 
initial condition F{0) = 0. The equation of motion for 
/(t,s) is given by dj{t,s) = [iE{t) + F{t)]f{t,s) i, 
where E{t) = Ei{t) — EQ{t). For this two-level atom 
and the systems in the following two examples, we have 
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FIG. 1. Fidelity dynamics of the two-level system with differ- 
ent 7 and Red solid lines for free dynamics, blue dashed 
lines for r — 2A, black dot-dashed lines for r = 3A, and 
green dotted lines for r = 6A. The other parameters are set 
as u = 0.2r, and A = 0.04rt. 



E(t) ^ Lj + c(i), where uj is the bare frequency for the 
system and c(t) is the control function. 

With the initial state of the system |i/'o) = (1/ V2) (|0> + 
|1)), our aim is to control the population of this chosen 
state. The exact stochastic transition amplitude is 



1+ / ds'z*e-'Io ds[Eiis)-Eo{s)] 



.(11) 



where E[{s) = Ei{s) — iF{s). Consequently in the rotat- 
ing picture of Hgys, we have 



l + Fl{t)+2FR{t)Fi{t) 



^0 



dsids2 



X a(si, S2)FR{si)FR{s2)K{si)/k{s2) 



(12) 



where FR{t) = e- ^^Fi^i^) ^ Fi{t) = cos[ j*^ dsFi{s)\ {Fr 
and Fj stand for the real and imaginary parts of F re- 
spectively) and K{t) = e^!od'>[J^i'>)+Pi(.'>)]. 

For simplicity, throughout the paper the non- 
Markovian environmental noise is described by an 
Ornstein-Uhlenbeck type correlation [l^ a{t,s) = 
i^g-7|t-s|^ where F is the environmental dissipation 
rate and I/7 characterizes the memory time of the en- 
vironment. The finite, but non-zero I/7 gives a non- 
Mar kovian process with the Markov limit when 7 — ^ cxi. 
The equation of motion for F{t) can be easily derived 
from the QSD equation: 

dtFlt) ^^ + [-^ + ^E{t)]F{t) + F\t). (13) 

Therefore, Eq. ([T^ may be controled by the functions 
E{t) (or c{t)) via Eq. 



When c{t) = 0, the two-level atom will be driven by 
the environment into the ground state after a period of 
time. The red lines in Fig[T] show that the fidelity ap- 
proaches 0.5 regardless of the memory times. Clearly, in 
the absence of control function, the smaller 7 typically 
gives rise to more significant fiuctuations, yet it does not 
effectively preserve the fidelity. However, when the con- 
trol pulses are applied, we see that the non-Markovian 
features can greatly enhance the efficiency of the rectan- 
gular pulses for controlling the state of interest. The blue 
dashed line in FiglTJa) indicates that the fidelity can be 
preserved for a long period of time. Particularly, when 
the period of the pulse r is chosen as small as twice the 
interaction duration time A in one period, the fidelity is 
fully preserved up to the time Ft = 13. Moreover, we 
show that the control pulses can work very well even if r 
is as long as the three times of A (the black dot-dashed 
line). But if the period of pulses is too long, the dy- 
namical decoupling becomes highly inefficient [20l |. For 
example, when r = 6 A, as shown by the green dotted 
lines in FigsllJa),[T](b) and[ljc), the control function c(i) 
becomes adversary of the fidelity control. In all those 
three cases, the fidelity decays faster than those of free 
evolution (i.e., c(t) = 0). With decreasing memory time 
of the environment plotted in FigsHJb) andHfc), the con- 
trol effect of the fidelity is gradually weakened even using 
the same rapid pulses as in Fig|IJa) and simultaneously 
raising the interaction intensity ^ for partial compensa- 
tion. When 7 = 2.0 and r ^ 3A, we see that the control 
becomes inefficient in the Markov limit (large 7). 



B. A qutrit dissipative model 
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FIG. 2. Fidelity dynamics of the qutrit system with different 
7. In (a) and (b), the approximate results (blue dashed lines) 
by Eq. (|14p is compared with the exact ones (red solid lines); 
In (c) and (d). Red solid lines for free dynamics, blue dashed 
lines for r — 2 A, black dot-dashed lines for r — 3 A, and 
green dotted lines for r = 6A. The other parameters are set 
as a; = l.OF, A = 0.04n, and k = V2. 
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For this model, we have H^ys = E{t){\2){2\ - |0)(0|) 
and L = k(|0)(1| + |1)(2|). In the basis |1), |0) and |2), 
the efltective Hamiltonian takes for the following form: 



off 



-iFi 







— inUz 

















E 


-iF2 



(14) 



where F2(t) and C/^(t) = /J ds'C/(t, s')z;, are the 

three coefficient functions in 0{t,z*) = Fi(t)|0)(l| + 
F2{t)\l){2\ + U,Xt)\Q){2\. One particularly interesting 
feature of this model is that in the PQ partitioning 
W 0, but D is diagonal, so it allows a simple solu- 
tion. Applying the noise-free approximation, we have 
O{t,z*)=Fi{t)\0){l\+F2{t)\l){2\ with 

dFi{t) T"fK^ 



dt 

dF2{t) 



2 



+ (-7 + iE)Fi + Fl - 
+ {~^ + iE)F2+Fi, 



F1F2, 



(15) 



dt 2 

and Fj{0) — 0, j — 1,2. The approximation validity in 
Eq. (ini) is testified by Figiaja) and FiglSIb) with two 
difi'erent 7's. We could also find the analytical expres- 
sion of the fidelity by PQ partitioning technique. When 
= (1/V3)(|0) + |1) + |2)), the compact formula for 
the fidelity is given by 

ds dsia{si,s) 

'Q Jo 

X [F*{t)B*{si)+F*{si)E*{si)][F*{t)B*{si) 

+ F^isi)E*{si)]+ [ ds [ ds' [ dsi [ ' ds2 
Ja Jo Jo Jo 

X [a(si, s)a{s2, s') -I- a{si, s')a{s2, s)] 



I + h{t) + F2{t)[' 



X F*{si)E*{si)B*{s2)Fi{s)E{s)B{s'] 



(16) 



where Fj{t) 



^!^dsF,(s)^ ■ ^ ^ E{t) = e-*/o 



&nAB{t) = E{t)F2{t)lFi{t). 

In Figsl^ljc) andHJd), we compare the control dynam- 
ics with the same environments as in FigslJJa) and[2jb) 
respectively. As in the first example, we also notice the 
larger environmental memory time (less 7) is very help- 
ful to relieve the requirement of the pulse frequency than 
the smaller one (bigger 7). If 7 = 0.5, r ^ 3A, the 
fidelity could be maintained above 0.9 even when Tt ap- 
proaches 40. In Figl^fd), 7 = 2.0, it is interesting to 
observe a damping enhancement of the fidelity especially 
when T ^ 3 A. This (as well as all the green dotted 
lines in Fig[l]) corresponds to a sort of Quantum Anti- 
Zeno Effect (AZE) [21]. AZE occurs when the evolution 
is repetitively interrupted by projecting the system onto 
some state as in a measurement process [22!|, by periodi- 
cally applying the pulses [2^ with long time intervals, or 
by coupling between the excited state and an auxiliary 
state [2J1. Here it is induced by unitary pulse sequences, 
which renormalize the frequency of the system with a 
proper period and mimic a repeated measurement pro- 
cess. 



C. A special (N + l)-level atom 




40 



FIG. 3. Fidelity dynamics of a 101-level system in a non- 
Mar kovian environment. Red solid line for free dynamics, 
blue dashed line for r = 2A, and black dot-dashed lines for 
T = 3A. The other parameters are set as a; = 0.2r, A = 
QMVt. 

For this special case, we consider a genuine multi-level 
atomic system represented by i?sys — T^^=o ^n\n) {n\ and 
L = En^i where Eo = ~E{t) and E^^o = E{t), 

and for this model 0{t, z*) — F{t)L. Therefore, 



/~E 



HcS — 



\ 



E~iF -iF 
~iF E-iF 



-iF 



-iF 



-iF 
-iF 

E-iF J 



(17) 



where Fit) satisfies dtF{t) ^ ^ + {-j + 2iE)F + NF^ 
and -F'(O) — 0. Due to the time symmetry of D, 
[D{t),D{s)\ =Q [t ^ s), the propagator in Eq. ^ can 
be exactly obtained as 



G(t,0) 



1 r ^ 

-|5][^(t)i?^(o + (A^-i)^(i)] \m 



Y,[mF''(t)-Emn){r 



(18) 



where F{t) = J'o ^(^^''^ For the initial state \i>o) = 

{\/\/N + 1) Y^n=o I"")' ^'g^'in the PQ partitioning allows 
to find an analytical expression for the fidelity function: 



Ht) = 



1 + N"^ Elicit) + 2NF^ [t)F[[t) + N 



X f f dsids2a{si,S2)Fg{si)F^{s2)K'(si)/K\s2)\, 
Jo Jo 

(19) 

where F"j{t) = cos[J* dsNIpis)] and K'it) = 
^iJ^ds[2Eis)+Ni,.{s)]^ It is easy to check that when iV = 1, 
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Eq. p9|) reduces to Eq. (IT2)) . It is interesting to consider 
the weak coupling limit when F ^ 0, then D becomes 
a diagonal matrix diag[i?, E, - ■ ■ ,E], thus Eq. ([7]) can be 
simplified as 



l + N^ + 2N + N^ 



"'0 



dsids2 



(20) 



We emphasize that our formalism works for an arbi- 
trary + 1 system. For example, we take N = 100 
shown in FiglJl Without the control function c{t), the 
fidelity will quickly decays to l/(iV + 1). Clearly, the effi- 
cient control can be made possible only when the system 
is far from Markov regime. As shown in this example, our 
pulse control scheme works well for the non-Markovian 
environment with 7 = 0.2. In fact, by setting ^ = 4.0, 
and r = 2A, we show that T{Tt = 40) can be maintained 
as high as 0.85. 



III. CONCLUSION 

We have derived an exact one-dimensional stochas- 
tic master equation based on the Feshbach PQ parti- 
tioning and the non-Markovian QSD equation and ap- 
ply it to the quantum control dynamics of three distinct 
model systems. The periodical rectangular pulses are 



used to protect the subspace of interest from leakage to 
the other part of the Hilbert space of the system and 
the environment. The control dynamics measured by 
time-dependent fidelity is realized by tuning the pulse fre- 
quency l/r with different memory time 1 of the non- 
Markovian environment. Our results have exemplified 
the simplicity and power of the exact one-dimensional 
stochastic master equation. Our approach has paved a 
way for studying the direct control dynamics of an ar- 
bitrary multi-level atomic system without invoking the 
exact or perturbative non-Markovian master equations. 
Our hybrid technique would be versatile enough to ac- 
commodate other types of quantum control if some neces- 
sary modifications on the non-Markovian QSD equations 
are made. It is also possible to consider the quantum 
control dynamics of multi-particle system ,251 based on 
the PQ partitioning. We leave these open questions for 
the future investigations. 
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